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Abstract. In this note we analyse an anisotropic, two-dimensional bootstrap per- 
colation model introduced by Gravner and Griffeath. We present upper and lower 
bounds on the finite-size effects. We discuss the similarities with the semi-oriented 
q ! model introduced by Duarte. 
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1 Introduction 



Bootstrap percolation (BP) models, sometimes also called k-core percolation or 
threshold growth models, are simple Cellular Automata with a percolation con- 
figuration chosen as a random initial condition. The deterministic dynamics is 
given by repeatedly applying a growth rule, and the question of interest is if the 
whole system will be occupied in the long run. They have been a topic of interest 
in fields as diverse as physics - for the study of magnetic models, metastability, 
rigidity, fluid flow in porous media and glassy dynamics -, mathematics, com- 
puter science, neural science and economics, see e.g. HI HH5H6M71 ISHU EH EH EH 

EZl EE EE US 1221 E3 M 123 M HI • 

Most of the models which have been considered in the literature are isotropic. 
In this paper we consider an anisotropic model in two dimensions which was in- 
troduced by Gravner and Griffeath ||T5l[T6ll . Their model is critical, meaning that 
there are infinite sets with an infinite complement which can fill the entire lattice, 
but no finite set can do so, and all sets with a finite complement will fill up the 
lattice. It has a trivial percolation threshold p c = 0. We determine the asymptotic 
behaviour of the finite-size effects, and find it to be different from the one pro- 
posed by Gravner and Griffeath. We discuss the similarity with the behaviour of 
the semi-oriented (Duarte) model which is also anisotropic, and the behaviour of 
which was established by Mountford. HT3ll2ll3l|29l|24l|25ll. 



2 Definition and some known properties of bootstrap 
percolation models 

We will consider BP models either on 1? or on finite subsets thereof. 

The dynamics is defined as follows: 

Once a site is occupied, it will remain so forever. 

For each empty site it is checked at each time step whether in a prescribed neigh- 
borhood x + TV of the site x at least 9 sites are occupied. (The parameter 9 is called 
the threshold value.) If so, at the next time step the site will be occupied. 
At the next time step this procedure is repeated. 

Thus this updating rule provides a deterministic, parallel dynamics. This dy- 
namics runs until no changes occur anywhere anymore. 

To describe the evolution of an occupied set AcZ 2 we define the discrete time 
step operator T: 

T(A) = AU {x : \(x+J\f)nA\ > 9}. (1) 
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Starting from some initial set A c Z 2 , we iterate: T(A n ) = A n+ i. We also define 
Aoo = T°°(A ) = U^Lo A n - We call a model critical if there exists some infinite A 
(but no finite one) with an infinite complement in Z 2 such that A^ = Z 2 , that is 
the dynamics eventually fills all of Z 2 . 

We choose for our random initial configuration a percolation configuration, 
where each site is independently occupied with probability p. 

The standard (isotropic) model is defined by taking as the neighborhood of a 
site its four nearest neighbours and 9 = 2. 
• 

Af = • • 

In that case p c = for the infinite lattice fl4| (so for each initial density the 
lattice will fill up). The finite-size effects at small p are such that if the size of a 
square is of order exp(Ci X-), with G\ > Cst for an explicit, computable constant 
Cst, the square tends to be occupied in the end, while if it is of order exp(C 2 x 
-), with C 2 < Cst, it tends to remain mostly empty BSHlZl with overwhelming 
probability. 

In the Gravner-Griffeath (GG) model, the neighborhood of a site consists of six 

sites, to wit, the closest two sites to the East, the closest two sites to the West, and 

the nearest neighbours in the North and South direction. The threshold value 9 

now is defined to be 3. 

• 

Af = • • • • 
• 

For later reference we also mention the Duarte (semi-oriented) BP model, in 
which the neighborhood of a site consists of three sites, namely its North, West 
and South neighbours, and 9 is chosen to be 2. 
• 

Af = • 
• 

This model can not grow in the Western direction. 

The GG model (similarly to the Duarte model) has easy and hard growth di- 
rections. Indeed, once an N by iV square is occupied, a single occupied site at 
distance 1 or 2 will fill the next line segment, in the East or West direction, while 
for growth in the North or South direction two sites with no more than 3 empty 
sites in between need to be occupied, which is much less probable (locally the 
probability is of order p 2 instead of order p). 

Two important notions in the theory are those of a "critical droplet" (some- 
times called nucleation droplet) and a volume being "internally spanned". 
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A critical droplet is an occupied configuration which will keep on growing 
with large probability. In the regular model, for example, an occupied square 
with side length C\ x -, C\ large, has this property, while in the Duarte model, 

an occupied line segment of length C\ x - x In -, C\ large, is known to have this 
property 

The second important notion is that of a volume being internally spanned. 
We will call a volume V 'internally spanned' if given some occupied set A c V, 
occupies all the sites in V. In that case the initially occupied set A - which 
typically will be a random set - spans the volume V. (In the Duarte model the 
notion of internally spanned is defined slightly differently.) 

Remark: A related quantity, for which a similar scaling behaviour can be de- 
rived, as a corollary of our results, is the first passage time at the origin, that is the 
time at which the origin is occupied for the first time, see e.g. I51H511. 

3 Main result and proof 

Our main result is that the volume size for which the Gravner-Griffeath model 
changes from typically not being internally spanned to typically being internally 
spanned scales as exp(0(^ x ln 2 (i))). 
Indeed we have: 

Theorem 1. Consider the Gravner-Griffeath model. 

A) If a square is larger than exp(Cx(i x ln 2 (i))), with C\ large enough, it is internally 
spanned with large probability ', for p small enough. 

B) If a square is smaller than exp(C 2 (- x In 2 (-))), with C 2 small enough, it is not in- 
ternally spanned with large probability, for p small enough. 

Proof: We will prove A) and B) separately. 

To prove A), we will first show that an occupied rectangle of size 2 by C\ x - x 
In i is a critical droplet. 

Indeed, let such rectangle (a double segment, long in the North-South direc- 
tion, width 2 in the East- West direction) be occupied, then we will first argue that 
with large probability it will grow in the Eastern and Western directions, until it 

fills up a rectangle of size 2(i) C by C\ x i x In i as long as C is smaller than C\. 
(When C is larger than 1, such a rectangle is much larger in the East- West than the 
North-South direction). 



4 



For this to happen, it is sufficient that in all the (-) c line segments on the left 
and on the right at least one site is occupied in the original configuration. But this 

happens with probability (1 — (1 — p) Cl p ln p) " ~ (1 — p Cl )°^p \ which is close 
to one when C is smaller than G\. 

Once the larger rectangle is occupied, with C chosen to be large enough (it 
should be at least larger than 2), it will grow in all directions, with large proba- 
bility. Indeed, if we take the rectangle to be Ni by N 2 , N± much larger than i and 

N 2 much larger than (-) 2 , then the probability that the rectangle does not stop 
growing in either horizontal or vertical direction is larger than 

nr=i(i-(i-pT 2+fc )(i-(i-pr i+fc )- 

Minus the logarithm of this probability is approximately 
Y^kLiO- ~~ p 2 ) N2+k + (1 — p) Nl+k which is small, showing that the corresponding 
probability is close to I. 

As the probability for a particular rectangle to be a critical droplet (thus the 
density of critical droplets) is larger thanp 2Cl p ln p = exp(— 2Ci-ln 2 (i)), the mini- 
mal size necessary for a square to contain with large probability a critical droplet 
is bounded from above by the inverse of this density. This critical droplet then 
will grow to fill up the entire square. Indeed, wherever in the square the droplet 
is, for small p there will always be enough space to grow sufficiently far in at least 
two perpendicular directions (remember that the size of the square is much larger 
than \) and once the occupied set is big enough, with large probability it will then 
fill up its complement in the square. 

This finishes the proof of part A). 

To prove part B), we essentially follow the analysis of Gravner and Griffeath 
lH~5l , and point out where it needs to be modified. (Gravner & Griffeath based 
this part of their analysis largely on the proof of the similar statement for ordinary 
bootstrap percolation as given by Aizenman & Lebowitz in [5].) 

We will first quote their definition and main claim: 

Definition: [We] call a rectangle RinZ 2 potentially internally spanned (PIS) if it 
is either a single site in II (p) or if (z) for every (integer) vertical i through R there 
exist two sites x, y e n(p) n R such that \\x — y\\oo < 4 and they are both at £°°- 
distance at most 4 from i, and (z'z) every horizontal line i through it has a site in 
II(p) fl R at £°°-distance at most 2 from t. 

Here II (p) denotes the set of initially occupied sites, when the occupation prob- 
ability is p. The property PIS is obviously weaker than that of being internally 
spanned. 
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Claim: Let L < M be positive integers. Assume that the origin is not eventu- 
ally occupied if the dynamics are restricted to [-L, L] 2 , but is eventually occupied 
if the dynamics are restricted to [-M, M] 2 . For every integer a e [4, L/4], there 
exists a PIS rectangle R included in [— M, M} 2 whose longest side is between a and 
4a. The proof to the claim is given in |fi~5l . 

From the definition of PIS it is then derived that the following must hold: 



P(k, I) := P(a fixed k x I rectangle is PIS) 

<(l_(l_ p 2 )fc) * (l _ (l _ p) y 

<min{(l-(l-p 2 ) fc y,(l-(l- P m 
:= min{/, II} 

(2) 



Now we deviate from [15J. We choose k = -gW, and I = C 2 - In-, then we have 



//= (l-(l-p)^ ln p)p 3l/2 

~ (1 _pC 2 ^^72 

~ exp 



p p 

if C*2 is small enough, and moreover 



p3/2 

<exp(-C 2 -ln 2 -) (3) 



/= (l_(l_p2)772)^-ln- 
~ (1 _ e -v^) C 4 ln ? 

p2p p 

= exp(-^ln 2 i) (4) 
2p p 

which shows that for this choice of k and I and Ci small enough 

P(M)<exp(-^ln 2 -). (5) 
2p p 
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This inequality contradicts the claim on the asymptotics of Gravner and Grif- 
feath. (We suspect that in their derivation, which they don't give completely but 
describe as a "straightfoward computation", the following went wrong: They 
seem to have used the approximation (1 — x) = exp(-xN) not only for small 

values of x, but also for larger x-values. For example, in an expression like I, with 

i 

the choice x = (1 —p 2 )p 3/2 , and N = C 2 x - x ln^, the approximation gives a 
wrong answer, which is of the form proposed in (l5|. For some further discussion 
on this point, and another derivation of part A) of the theorem, see (HI.) 

Thus if the size of a square is less than the inverse of this probability, it will not 
be PIS with large probability This is the main ingredient to prove B). For further 
details see again ||T5l . 



4 Discussion and conclusions 

We have found that the asymptotics of the finite-size effects is different from the 



one proposed in IH5I1 . In fact, the behaviour of the GG model is quite similar to 
that of the Duarte model. 

It has the same asymptotics, and whereas the Duarte model has a single line 
segment of length C\ x - x In -, with C\ large enough, as its critical droplet (first 
observed by Schonmann and communicated to the authors of (2l HI), in the GG 
model a double segment of the same length plays the same role. Therefore the 
difference between the two models suggested in fl5f (one model understandable 
via its critical droplets, the other one via its growth mechanism) is in our opin- 
ion spurious. This observation leads us to expect that in considerable generality 
critical anisotropic BP models with slow and fast directions will have the same 
asymptotics for their finite-size effects. 

Open questions: 

It would be interesting to see if it is possible to identify the exact asymptotic con- 
stant (in other words, can C\ and C 2 be chosen to be identical, as happens in the 
isotropic model ffTTI ?). 

Also it seems a natural question to see what the behaviour of anisotropic mod- 
els in higher dimensions might be. 

However, at present we have no results pertaining to these questions. 

Acknowledgements: We thank Christof Kiilske for a critical reading of the 
manuscript. 
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